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Abstract 

We consider here the dynamics of some homogeneous and isotropic cosmolog- 
ical models with N interacting classical scalar fields non-minimally coupled 
to the spacetime curvature, as an attempt to generalize some recent results 
obtained for one and two scalar fields. We show that a Lyapunov function can 
be constructed under certain conditions for a large class of models, suggesting 
that chaotic behavior is ruled out for them. Typical solutions tend generically 
to the empty de Sitter (or Minkowski) fixed points, and the previous asymp- 
totic results obtained for the one field model remain valid. In particular, we 
confirm that, for large times and a vanishing cosmological constant, even in 
the presence of the extra scalar fields, the universe tends to an infinite diluted 

matter dominated era. 
95.10.Fh, 98.80.Cq, 98.80.Bp 

Typeset using REVTgX 



*e-mails: lbrenig@ulb.ac.be, annibal@fis.unb.br, egunzig@ulb.ac.be, marciano@fis.unb.br, 
asaa@ime.unicamp.br 



1 



I. INTRODUCTION 



We have considered recently the homogeneous and isotropic solutions of the cosmological 
model described by the action [1]: 

S = \J d A x^~g {-(1 - £/# 2 )^ + gTd^W - 2V^)} , (1) 

where ip is the non-minimally coupled scalar field, k = 8nG, G is the Newtonian constant, 
R the scalar curvature of spacetime, and the self-interaction potential V{jp) has the form: 

V W = - ^ - - ? , (2) 

with a = m the mass of the scalar field, f2 is an arbitrary constant and — ^ is the usual 
cosmological constant A. The nonminimal coupling term in (1) is, of course, (,Rip 2 where 
£ is an arbitrary constant. The results previously obtained (see [1] for the motivations and 
references) were worked out for the case £ = 1/6, the so called conformally coupled case, 
and they point toward some novel and very interesting dynamical behavior: superinflation 
regimes, a possible avoidance of big-bang and big-crunch singularities through classical birth 
of the universe from empty Minkowski space, spontaneous entry into and exit from infla- 
tion, and a cosmological history suitable for describing quintessence in principle. Through 
exhaustive numerical simulations and with some semi-analytical tools, the 3-dimensional 
phase space of the model has been constructed. The existence of a Lyapunov function for 
the fixed points was crucial for the study of the asymptotic behavior [2] and for precluding 
the appearance of any chaotic regime, confirming and shedding some light on some other 
previous results in this line [3]. 

The robustness of its predictions must be an essential feature of any realistic cosmological 
model. The study of some generalizations of the model (1) is, therefore, mandatory. In [4], 
we consider the implications of the inclusion of a second interacting massless scalar field 
in the dynamics of homogeneous and isotropic solutions. The corresponding phase space 
becomes now 5-dimensional, and richer structures could appear. In spite of this, we show 
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that for a class of physically reasonable quartic interaction potentials, the neighborhoods of 
the relevant fixed points are unaltered, and, in particular, the asymptotic regimes obtained 
for the 1-field case [2] are preserved in the presence of an extra massless field. The relaxing 
of isotropy was the next performed robustness test of the model. In fact, it appeared as 
a corollary of a much more general result. In [5], we studied the singularities in the time- 
evolution of homogeneous and anisotropic solutions of cosmological models described by the 
action: 



with general F(ip) and V(ip). Such a model has been recently considered also in [6]. We 
showed that anisotropic solutions generically evolve toward a spacetime singularity that 
corresponds to the hypersurface F(ip) = in the phase space. Such singularity is harm- 
less for isotropic solutions. A second, and different, type of singularity corresponds to the 
hypersurface Fiijp) = 0, where 



and no solution can avoid it. Starobinsky [7] was the first to identify the singularity cor- 
responding to the hypersurfaces F{jp) = 0, for the case of conformally coupled anisotropic 
solutions. Futamase and co-workers [8] identified both singularities in the context of chaotic 
inflation in F(ip) = 1 — theories (See also [9]). For this type of coupling, the first singu- 
larity is always present for £ > and the second one for < £ < 1/6. The results of [5] are, 
however, more general since the case of general F(ip) is treated and all conclusions are based 
on the analysis of true geometrical invariants. The phase space for the model (3) is also 
5-dimensional, and many new structures appear. The appearance of the singularity of the 
first type implies the instability of the anisotropic solutions even for the conformally coupled 
case, for instance. The presence of any amount of anisotropy, no matter how small, makes 
the model unstable, challenging its validity as a realistic cosmological model. The asymp- 
totic behavior in the neighborhood of the fixed points far from the hypersurface F(ip) = 
is, however, preserved. 




(3) 



F 1 W=FW + 3 -(F'W) 2 , 



(4) 
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We present here a new robustness test for the homogeneous and isotropic solutions for 
models of the type (1). We consider the case of several interacting scalar fields conformally 
coupled to spacetime curvature. Exact solutions for this case are much more difficult to find, 
in particular we could not identify the heteroclinics and homoclinics as it was done in [1]. 
We could, however, get a strong result about the stability of fixed points suggesting that 
chaotic regimes are ruled out. This is proved for a class of interaction potentials, through 
the construction of an explicit Lyapunov function. In spite of the unstable character of 
anisotropics for models of the type (1), the results reported here are expected to be valid in 
the neighborhood of the Minkowski fixed point, even for the anisotropic case. 

II. THE LYAPUNOV FUNCTION 

Let us consider the action for iV classical interacting scalar fields ipi, . . . , ipN conformally 
coupled to the spacetime curvature: 



S=\j d'x^-g {-[1 - + ... + i> 2 N )\* + £ jT^M^i - 2V(^, . . . , <M j 



(5) 



where the potential V is a polynomial in the N fields variables ipi up to fourth degree. We 
consider homogeneous and isotropic solutions corresponding to Robertson- Walker metrics 
with flat spatial section 

ds 2 = dr 2 -a 2 (r)(dx 2 + dy 2 + dz 2 ). (6) 

Variations of the action (1) with respect to the scalar fields fa yield N coupled Klein- 
Gordon equations 

• 1 dV 

^ + 3HiPi - -RiPi + — = 0, (7) 
6 dipi 

where H — - is the Hubble function and the dot denotes time derivative. The variation 

a 

of S with respect to the metric leads to the Einstein equations which can be cast, in an 
analogous way to the one field case, in the form 



4 



K N K N N 

Y,tf + *V- 3H 2 + -H 2 ^l + #E^ = 0. 

i=l " i=l i=l 



(8) 



i.e. the energy constraint, and 



/ N dV\ 

R = -Q(H + 2H 2 ) = k \-W + ^ ^_ J , 



(9) 



the trace equation. The system of ordinary differential equations (7) and (9) is defined 
in a (2iV + f )-dimensional phase space. The solutions are, however, confined on the 2N- 
dimensional zero energy hypersurface (8). 

For physical interpretation reasons (the same ones of [1]), we assume that the potential 
V(ipi, . . . ,iPn) contains a mass contribution in the form of a quadratic term for each field 
■ipi, and a homogeneous quartic function f^ipi, . . . , ip N ) describing the interactions of the N 
scalar fields (including possible self-interactions), and a cosmological constant 



3 N 9lo 

v(^i, . . . , i/> N ) = - E + fa(*h.> • • • , M — J 



K 



i=l 



K 



(10) 



where cej = |mf , being the mass of field ipi. 



We now proceed to show the existence of a Lyapunov function for the dynamical system 
(7)-(9) under some well defined conditions on the potential (10). We focus on the sub-system 
of Klein-Gordon equations (7). By multiplying each equation by ipi and summing on the N 
fields we get 



d 



N ih 



dr 2 



N 



N 
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(ii) 



i=l 



Assuming that —^Ripi derives from a potential 
we would obtain: 



(12) 



d_ 

d^ 



N il 2 

Y,f + v + u 

1=1 z 



N 



(13) 



i=i 



If the function between square brackets has a minimum at a given fixed point, it will be 
a candidate for a Lyapunov function for H > 0. Equations (7)-(9) have, obviously, many 



fixed points, but we are concerned only with de Sitter (cu < 0) or Minkowski (cu = 0) 



fixed points, for which ^ = — and H = ±y — "iuj / k. We are restricted, therefore, to 
uj < (non-negative cosmological constants, as in [1,2]). A closer analysis of the energy 
constraints (8) reveals that the system can cross the 2iV-hypersurface H = only on a 
27V" — 1 submanifold of the 2N + 1 original phase space. In a neighborhood of the fixed point 



(H = y — 3u/k, ipi = 0) where V is non-negative, the system can reach H = only at the 
point where V vanishes, otherwise the system cannot reach H = 0. Hence, we will have for 
H > 



d 

dr 



2 



tr + V + U 



< (14) 



in such a neighborhood and, provided that V + U has a minimum at the origin, 

N i) 2 

L = Y,f + V + U (15) 

i=l 

is a Lyapunov function, ensuring, thereby, the stability of the fixed point (H = 



^-3u/K,x/ji = 0). 

The relevant question here is the validity of the hypothesis (12). It is valid, for instance, 
if all the fields have the same mass, a.\ — a-i — . . . — a^- This is a strong constraint that, 
as we will see below, can be somehow relaxed. From equation (9), with the potential given 
by (10), we have 



Condition (12) requires 



■^Ripi = ipif2 <Xji/>j ~ —ipi (16) 
o j=1 K 



d(R^) _ d(R4> k ) 



implying, from Eq. (9), that 



otk^k = aii>ii>k (18) 



for all couples (i, k) with i ^ k. The equal masses case is, obviously, a particular solution of 
(18). For this case, the potential U is readily shown to be 
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Ufa, v> 2 , ■ ■ ■ , M = u + j\i>\ 4 - — |vf (19) 

where C/o is an arbitrary constant and | | stands to the usual Euclidean norm 

|Vf = ^ + ... + V^. (20) 

In this case, the Lyapunov function is given explicitly by 

L = ^\ 2 + -(a- c)H 2 + f 4 fa, ...,4> N )-^ + U + ^ (21) 

for arbitrary homogeneous functions of degree four / 4 . In summary, the function L has a 
minimum on the fixed point and there is a neighborhood (the attraction basin) with H > 
where L has non-positive time derivative. As the system can reach the hypersurface H = 
only in the fixed point, this provides a stability proof for the fixed point ipi — i>i — 0, 
H = +\j3\u\/n. An analysis of the time- reversed system reveals that the symmetric fixed 



point (with H = —<\J3\u\/k) must be repulsive. The trajectories starting in the vicinity of 
the last point will leave this region and some of them will cross the hypersurface H = 0. 
Some of the crossing trajectories can be eventually trapped in the half-space H > and 
will tend asymptotically toward the stable fixed point. This phase portrait suggest a regular 
behavior, and it is very reminiscent of the behavior exhibited by the solutions of the one-field 
model, where chaotic regimes were ruled out. 

As it was already noted, this reasoning breaks down for N scalar fields with distinct 
masses. Nevertheless, we now prove that, under rather weak conditions relating the masses 
and the parameters for some potentials of the form (10), another Lyapunov function exists 
and ensures the stability of the fixed point. We will consider here the special case where the 
homogeneous function / 4 has the form (10) 

N O N ( ft \ 

1=1 4 k,i=i(k&) V / 
with Pi k = (3ki- The Klein-Gordon equations (7) become 

^ + ZHiPi + - u)i/>i + {a z - + F i = (23) 

K 



where 



N N N 

j=iO¥0 j=Hj¥=i) j=i(j¥=i) 



(24) 



Let us now rescale the field variables by a constant factor ipi = jiipi. Klein-Gordon equations 



are now 



i>i + 3Hip l + -(<%- u)A + rMi - n^f + F t = o 

K 



(25) 



with 



N 

1 % 



i ' 



A? 



AT 



(26) 



We again multiply both sides of (25) by ibi and sum over i leading to 

■„2 



d 



N ib 



N _ • N -2 

i=l i=l 



(27) 



where 



A? 



i=i 



(28) 



N 



Let us determine the conditions under which the term ^Fjipj is a total derivative 

i=i 

This requires that Fj derives from a potential V 2 , that is 

5£ = 55 (29) 

Conditions (29) lead to the following algebraic relations between the masses, the coefficients 
of the potential V ad the scaling factor 7^ 



7 t 2 («i + A fe ) = 7fc(a fc + /3 ifc ), 
and = 0. The conditions (30) provides sign conditions 



(30) 



Oik + Pik 



> 



(31) 
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(34) 

= 1 



and homogeneous algebraic linear equations for the 7? whose compatibility conditions relate 
directly the masses with the interaction coupling constant (3^. For N = 2 (two scalar 
fields) the condition (30) are rather weak 

^±|H > 0. (32) 

«2 + Pl2 

For N = 3, these conditions provide the inequalities (31) and one strict equality : 

(ai + /3i 2 )(a 2 + ffaXas + &1) = 1 ^ 

It is important to remind that = f3ji. For N = 4, there are four independent conditions 

(Ol + /3r2)(c*2 + /^23)(Q3 + /?3l) _ j 

(«i + ^13) ("2 + #21) (a 3 + ^32) 

(q 2 + ff 23 ) (Q3 + ^34) (g4 + ^42) = 1 
("2 + Am) ("3 + ^32) («4 + A13) 

(ai + /3l 2 )(«2 + $24) (04 + flu) 

(«1 + /9l 4 )(«2 + /9l2)(«4 + $24 ) 

(ai + Pu)(a 3 + /3 3 i)(« 4 + gjg) _ 1 
(3 13 )(a 3 + (3 M )(a A + (3 41 ) 

There are four conditions and eight parameters ai,flij. More generally, the total number of 

parameters ctj, (3^ is ^ and the total number of conditions (strict equalities) is N ( N ~ 1 )( N - 2 ) % 

We, thus, can determine the maximum number of interacting scalar fields above which the 

above conditions are overdetermined 

tP N(N-l)(N-2) 

2 6 1 ' 

which yields the solution N > 5. Finally, the explicit expression of the Lyapunov function 
in the original variables ipi is 

E wf^^W (36) 



with the conditions (30) over the parameters ai,/3ij, and •yf. 

Again, this function vanishes at the fixed points and on the H-axis. It is non-negative 
definite in a finite domain of the phase space around that axis and its time derivative is 
given by 

§ = -3ffE4# (37) 

i=l H 

which is non-positive in the half-space H > 0. This proves the stability of the fixed point 
i>i — i>i — 0, H — All the solutions tend to this point or go to infinity Furthermore, 

both Lyapunov functions may be global, i.e. 

lim L = +oo (38) 

when the parameters satisfy some inequalities. In that case, the Lyapunov stability theory 
shows that solutions are bounded and accumulate on one of the fixed points, depending on 
the initial conditions. 

III. CONCLUSION 

We have shown under which conditions some N non-minimally scalar fields homoge- 
neous and isotropic cosmological models admit a Lyapunov function for their de Sitter (or 
Minkowski) fixed point. The physical interpretation of such conditions is still unclear. 
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